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ABSTRACT 

We present a spectroscopic deprojection analysis of a sample of ten relaxed galaxy clusters. We use 
an empirical F-test derived from a set of Markov Chain Monte Carlo simulations to determine if the 
core plasma in each cluster could contain multiple phases. We derive non-parametric baryon density 
and temperature profiles, and use these to construct total gravitating mass profiles. We compare these 
profiles with the standard halo parameterizations. We find central density slopes roughly consistent 
with the predictions of ACDM: —1 < (ilog(/3)/(ilog(r) < —2. We constrain the core size of each cluster 
and, using the results of cosmological simulations as a calibrator, place an upper limit of ~ 0.1 cm^ 
g^^ = 0.2 b(GeV/c^)^^ (99% confidence) on the dark matter particle self-interaction cross section. 

Subject headings: X-rays : galaxies: clusters — cosmology : dark matter 



1. INTRODUCTION 

The ACDM theory of cosmological structure forma- 
tion - a cold dark matter (CDM) universe evolving un- 
der gravitation with a cosmological constant (A) - has 
been very successful in characterizing large scale struc- 
ture on scales greater than ~ 0.1 Mpc (Mohr et al. 200S 
i Lahav et al. 2001; iSperjgel et al. 2003; Pcrcival et aT| 
120031 . However, ACDM makes three predictions which 
appear to conflict with observations. The first of these 
discrepancies is the "satellite problem" - the Milky 
Way contains at least an order of magnitude fewer 
satellites than simulations predict f or a da rk matter 
halo of its size (Kauffma n. White, fc Guideroml ITool 
IMoore et alllTQQQa^ .Klvpin et alJ 119991 The second is 
the "disk problem," - disk galaxies produced in simu- 
lations contain too little mass and angular momentum 
(|Navarro & Steinmetz 2000). The third is the "core 
problem" - galaxy-sized dark matter halos produced in 
simulations tend to have significantly steeper logarithmic 
density gradients in their cores than do observed struc- 
tures. 

The density profile of a dark matter halo formed 
through the hierarchical assembly of smaller structures 
is usually parameterized as a pair of power laws with a 
transition scale radius Ts fe.g. lJing fc Sutol l)2002D ). 

^ = \ (1) 

Po (r-/r«)"(l-hr/r-,)T-« ' 

with a = 1 ijNavarro. Frenk fc WhiT^ l|199fl'.'1997^: here- 
after t he "NFW" profile) or a = 1.5 (Moore et al. 
(|l,999b); hereafter the "Moore" profile). While there is 
general agreement that 7 = 3, rotation curves of low 
surfa ce brightness l|Burkerti 119951 iMcGaueh fc de Blo^ 
19981) and spiral galaxies ijSahicci fc BiirkertI '200(T: 
Gentile et alJl200Cl[) and the funda mental plane of ellipti- 
cal galaxies l|Borriello et alJl2003(l suggest that core den- 
sity profiles are significantly fiatter than ACDM, with 
> a > 1. 

^ Center for Space Research, Massachusetts Institute of 
Technology, Cambridge, MA 02139; jsaaspace.mit.edu, 
mwb@space.mit . edu 



Reports of a core problem on galaxy cluster 
scales have been somewhat more controversial. 
Some lensing studies have fouri d flat c ore profiles 
iTvson. Kochanski fc DeirAntoniol 119981 iSmail et all 
l2000() wh ile others, makiiig use of similar data sets, 
have not iBroadhurst et al.ll20 0G^ iShaoir o fc 111611200011 
Most recentlv. lSand. Treu fc E llis (20021 and lSand et alJ 
((2,004) have found shallow (p ~ ^-o.s-j ^ores in several 
clusters using a combination of radial and tangential 
strong lensing arcs in tandem with stellar velocity 
dispersion measurements. 

The severity and cause of these discrepancies 
have been a subject of intense debate for nearly 
a decade. While observational limitatio n s ma y 
play a role iSwa,ters. Madore fc TrewheO I2001D. 
they are unlikely to be the sole cause i n the 
majority of cases JPalcanton fc BernsteiiJ 120001: 
Igwaters. Madore fc Trewhella' I200H iSwaters et^ 
12003; Bolatto et al. 2003). Various explanations have 
been offered which can be ascribed to astrophysics or 
to particle physics. Most ACDM simulations evolve 
collisionless dark matter particles and baryons under 
the influence of gravitation only; final baryon distribu- 
tions are usually determined according to prescription 
afterward. Thus, the disagreement between simulation 
and observation is either caused by neglected baryonic 
processes (astrophysics) which could act to fiatten 
the baryon-dominated cores, or could be the result 
of neglected processes involving dark matter (particle 
physics). 

It is certainly plausible that many of the problems 
on scales of ^ 10 kpc are the result of neglected 
baryon physics, if only because this is the scale at 
which star formation and the associated processes dom- 
inate the energy budget in galaxies. It seems less 
likely that discepa ncies ove r the range 100-800 kpc, 
scales probed in the lSand et al. ( 2004) study, could have 
the same provenance. Thus if these measurements are 
correct they may require a significant revision of our 
understanding of structure formation and the behav- 
ior of dark matter. The literature already provides 
a plethora of alternative dark matter theories: self- 
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interacting dark m atter llSpergel fc Stein hardt) l200Clj) . 
warm dark matter ll Hogan fc DalcantonI 120001) . annihi - 
lating dark matter ijKa plinghat. Knox fc Turnei^ '2000*) . 
scala r field dark matter (Goodman 2000; Hu fc Peebles 
]2Q0(f), decaying dark matter ( Gen 2001), and mirror mat- 
ter ( Mohapatra. Nussinov fc Teplitz 2002.; Foot 2004) . 
to name several. Of these, self-interacting dark matter 
(SIDM) ISpergel fc Steinhardt, (,2000|) has probably gar- 
nered the most attention, since with only one or two 
scattering events per particle over the age of the uni- 
verse the cores of galaxies could be sufhcient l y soft ened to 
match observed galaxy profiles IPave et aP l|2001f) (here- 
after DSSW). SIDM is also theoretically attractive in 
that i t would produce soft ened cores on galaxy cluster 
scales ^^shida et al. (2000), potentially providing an ex- 
planation for the Sand et al. (2004) results. 

The 10^ K X-ray emitting cluster plasma provides 
an important alternative probe to gravitational lens- 
ing reconstructions of cluster mass profiles. Lensing 
studies, which measure projected mass, may preferen- 
tially select clust ers with significant substructure along 
the line-of-sight l)Cohn fc Kochan'eS t2003). In addition, 
most lensing-deriv ed masses are mo del- dependent re- 
constructions (e.g. iSand et alJ l|2004H ). On the other 
hand, while X-ray mass determinations can utilize non- 
parametric spectroscopic deprojection, they are valid 
only if the hot plasma is in hydrostatic equilibrium. 
Those clusters for which this is most likely true usu- 
ally display "cooling flow" emission which can alt er the 
derived profile depending on how it is modelled f AUenl 
pjj98) ; Arabadiis, Bautz fc Garmirc (2002); Smith et aTj 
|2004 : lArabadiis. Bautz fc Arabadiis! |200^, hereafter 
ABA). Since these methods employ completely indepen- 
dent assumptions and potentially suffer from different 
systematic uncertainties, they constitute a powerful two- 
pronged attack on the core problem. 

It has been argued ijMarkevitch et alJl200,3t) that the 
complex cluster surface brightness distributions revealed 
by Chandra throw the hydrostatic hypothesis into doubt. 
We take a more optimistic attitude, noting that we select 
clusters which are largely devoid of significant substruc- 
ture. In addition, our method includes a diagnostic for 
clusters whose cores are not in hy d rostat ic equilibrium 
(see ©• We note that iKav etldl l|2004D find that ne- 
glecting a kinematic pressure term from a hydrostatic 
analysis leads to errors in the determination of cluster 
masses of only ^ 10% for those whose cores contain cool 
plasma. 

In this paper, we use X-ray imaging spectroscopy from 
the Chandra X-ray Observatory to constrain cluster mass 
profiles under the assumption of hydrostatic equilibrium. 
We briefly outline our mass profile extraction technique 
in fjS] We then examine the core mass profiles of a sample 
of relaxed clusters and compare with the predictions of 
ACDM in ^ Finally we use these profiles to constrain 
the self-interaction cross section of dark matter particles 



2. EXTRACTION OF MASS PROFILES 



The spectroscopic deprojection and mass recon- 
struction technique we employ has been ful l y de- 
scribed elsewhere fArabadiis. Baut z fc Garmird 120021: 
lArabadiis. Bautz fc Arabadiis ,2004) so we only briefly 
summarize them here. Starting with a Chandra ACIS-S 



or ACIS-I observation of a cluster, we process the data 
in the usual way^, removing point sources and periods of 
high background. We divide the events into concentric 
annuli centered on the emission peak such that there are 
at least 1000-1500 source photons per annulus. We con- 
struct a model consisting of a set of concentric spherical 
shells whose inner and outer radii and Tq correspond 
to the angular radii of the annuli in the data. Each shell 
is modelled as a thermal plasma at redshift z using the 
MEKAL model (lMewe^_GTo nenschild fc van den OordJ 
1986; iMe^^JLemei^^ van den Oord 1986; Kaastr^ 
199i iLiedahl. Ost erheld fc Goldstein 1995) within 
XSPEC l)Arnaiid| Il996i) . Each plasma component is 
characterized by a temperature T and normalization K . 
(The metal abundances are set to 0.3 solar.) The nor- 
malization is related to the plasma density p = n mp nu 
through the shell geometry. 



K-'^ = 1Q^^[ATTDA{1 + Z)f I UeUH 

J n 



dr (2) 



where Da is the angular diameter distance of the cluster. 

We test each cluster for the presence of a second 
(cooler) cospatial emission component in the core using 
an empirical F-test applied to a suite of 1000 Markov 
Chain Monte Carlo (MCMC) simulations per cluster. We 
use the MCMC significance S - defined as one minus the 
fraction of simulations which show an improvement in 
the model fitting as large as that shown by the data sim- 
ply by chance - to decide whether a second core emission 
component is required, adopting a threshold of 0.99 for 
its inclusion (ABA). 

We select clusters from the Chandra archive which dis- 
play (i) a single peak in their X-ray emission and circular 
isophotes, and (ii) no obvious signs of merging activity 
or other asymmetry in their X-ray surface brightness dis- 
tribution. Since these criteria are somewhat subjective, 
we further screen clusters in the sample according to the 
shape of their core mass profiles. Specifically, when de- 
riving ensemble averages, we use only those clusters for 
which a power law is an acceptable fit {Xr < 1-3) to the 
core profile (the innermost 5 points in the profile, usu- 
ally r < 100 kpc). We shall henceforth refer to these as 
the "relaxed core cluster" (RCC) subset. This prevents 
contamination of our sample from clusters which appear 
relaxed but whose irregular core mass profiles may belie 
the hydrostatic hypothesis. Seven of the 10 clusters in 
our original sample (ABA) are RCCs, though we include 
the remaining three objects for comparison in parts of 
this analysis. TableHlists the 10 clusters in the sample, 
their redshift, the treatment adopted for their core emis- 
sion, the MCMC significance of the second component, 
and whether or not the cluster core is relaxed. 

The model parameters, Ti, K {i = 1, N), and a Galac- 
tic absorption column A^hi are fit for simultaneously, and 
2D confidence limits for each Ti, Ki (alternatively Ti, pi) 
pair are calculated. The runs of temperature and den- 
sity for each cluster are then used to calculate the total 
enclosed gravitating mass of the c luster Mr using the 
hydrostatic equation l|Sarazinlll988D : 

M. = 3.535xl0i"M, ( ^ \ ( ' \ ( ^^"^^^ ^^^^^ 



IkeV 



1 kpc 



d log r d log r 
(3) 



See |http: / /asc/harvard.edu / ciao / threads /all, html | 
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Statistical fluctuations in the temperature measurements 
can result in an unphysical mass determination, i.e. 
Mr < or dM/dr < 0, over a limited range in r. To 
remedy this situation we impose the constraints Mr > 
and dM/dr > on the mass profile and allow the tem- 
perature and density of each shell to vary slightly in order 
to satisfy them. The degree to which these constraints 
are enforced is a free parameter (i.e. a weighting factor 
in a penalty function); if the variation in T, or less likely 
p, required to satisfy the constraints is excessive (e.g. 
1ST I ot ^ 1) we conclude that one of our input assump- 
tions, and therefore the mass profile, is suspect (ABA). 

3. CORE PROFILES 

In Figurenwe show mass profiles (top) and cumulative 
power law slopes a (bottom) for all 10 sample clusters, 
with Uniphase models at left and two-phase models at 
right. The exception to this is Hydra A, where we show 
only a Uniphase model. Our attempts at fitting a sec- 
ond core component in this cluster resulted in its best-fit 
temperature m atching that of the original component, 
consistent with lDavid et al. (200^. For each cluster, the 
model with the higher MCMC F-test significance (ABA; 
see also Table is indicated by the 'preferred' label. 

The cumulative power law slope, plotted as filled 
squares, represents a power law fit to the profile as a 
function of the outermost point used in the fitting. The 
innermost value of a is duplicated at the minimum r 
value with an open square for clarity (except for Hydra 
A and A1795, whose core profiles are unphysical even 
with a significant weighting of the mass constraints term 
in the penalty function). ACDM predicts that a should 
asymptotically approach —1 (the NFW profile) or —1.5 
(the Moore profile) as r decreases. 

3.1. Uniphase versus two-phase 

In general two-phase models have steeper central mass 
profiles than do uniphase models. This is so because 
the two-phase models allow the temperature of the hot 
component to remain roughly constant all the way in to 
spatial resolution limit. This augments the mass at each 
radius (see EquationEl, making Mr flatter in the center. 
Since Mr ~ j-"'^^, this implies a steeper density profile. 

3.2. Comparison with previous studies 

The profiles of power law slope show variation among 
the sample clusters, and in some cases between dif- 
ferent plasma emission models (Figure 0. The core 
reg ion of two of the clust e rs in our sample - Hydra 
A llMcNamara et alJ l2000t iDavid et~all l200l . A1795 
ijMarkevitch et alJl200H lEttori et al.ll2002tl . are known 
to be out of hydrostatic equilibrium. Figure ^ shows 
that the inferred mass profiles for these clusters are in- 
deed irregular. Three of the most relaxed clusters in our 
sample - A2029, A1689 and MS2137 - have been mod- 
elled in detail in other studies, allowing us to compare 
mass profiles. 

3.2.1. A2029 

The inferred core profile of A2029 is sensitive to the 
assumed plasma model (Figure^). The uniphase model 
has a nearly flat core, whereas the two-phase core, which 
is preferred > 99% significance, is consistent with the 



Moore profile. Confidence contours in a and the nor- 
malization Mq (where Mq = A/r(lMpc)) are shown in 
the upper left panels of Figure [3 The preferred two- 
phase model shows a core density that is even steeper 
th an the Moore profile. T hese results differ from those 
of iLewis. Stocke fc Buotel l|20021 . who find no evidence 
for multiphase core plasma, based on a lack of signifi- 
cant reduction in the value of their spectral fits, and 
Lewis. Buotc & Stocke (2003) (hereafter LBS), who find 
that NFW provides a good fit to their mass profile. These 
discrepancies may be due in part to the specifics of the 
deprojection method used in each case. Our technique 
requires more photons per annulus than the "onion peel- 
ing" approach, essentially because it is non-parametric 
and global, to adequately constrain the temperature at 
each radius. We must therefore use larger regions for the 
extraction of spectra. The result is that we resolve tem- 
perature structure down to about 15 kpc, whereas LBS 
resolve it to ~ 3 kpc. In their model the plasma tempera- 
ture drops a factor of two over this range, suggesting that 
perhaps the strong two-phase signature we see is due to 
a unresolved core temperature structure. We find strong 
evidence, however, for a second plasma component out 
to 34 kpc, well beyond our (and their) resolution limit, 
whereas LBS find none, so this cannot be the sole cause 
of the disparity. We note that Clarke et al. (2004) also 
find evidence for a second (cooler) emission component 
within r < 40 kpc. 

In addition to spatial resolution, a fundamental dif- 
ference between these two studies is the metric used to 
gauge the significance of the second component. LBS 
conclude that there is no evidence for multiphase plasma 
based on th e lack of significant reduction in the value 
of their fit l|Lewis. Stocke fc Buotdl2002D. Our metho d 
adopts the recommendation of iProtassov et all 1)2002(1 . 
who have studied statistical tests used to compare pairs 
of models when one lies on a parameter-space boundary 
of the other (in our case, the normalization of the sec- 
ond component set to zero). The have found that the 
standard statistical tests can underestimate the signal 
in noisy data. In such cases one should use a statisti- 
cal test performed on an empirical distribution derived 
from a large sample of data realizations which follow 
the complete paramet er likelihood function. Following 
IProtassov et alJ l|2002j) . our value for the significance of 
the second component is derived from the F distribution 
of a large MCMC sample (see We believe therefore 
that our value for the the core slope of A2029 is more 
reliable. 

3.2.2. A1689 

lAndersson k Madeiskl l|2004D used XMM-Newton 
data and the deprojection scheme of ABC to derive a 
core slope of —1.27, which is consistent with our results. 
It should be noted that they find evidence for an asym- 
metry in the temperature structure of the plasma, and 
suggest that this could be due to a merger event. Our 
mass profile shows no evidence of a merger, although our 
model does show a temperature jump at this distance 
(ABA, Figure 5). It is not clear that this affects the 
value obtained for the core slope, however. 

3.2.3. MS2137 
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Our profile for MS2137 is radicall y different from the 
strong gravitational lensing results of lSand. Treu fc Elhsl 
l)2002j) . These authors use radial and tangential arcs in 
conjunction with the stellar velocity dispersion profile of 
the central galaxy to constrain the central density pro- 
file, finding a > —0.9 at the 99% confidence limit. Our 
analysis yields a = —1.6 ± 0.2 (Icr), with a < —1.2 at 
the 99% CL. 

This discrepancy cannot be due to the differ- 
ence between deriving a dark matter profile profile 
ijSand. Treu fc Ellis 2002.) and a total matter profile (this 
work). Baryons dominate the mass profile only up to 
10 kpc l|Sand et alJl2004) : the innermost of our mass 
measurements is a t r ~ 20 kpc. In addition, the 
ISand. Treu fc Ellis! l|200^ analysis used a fixed NEW 
scale radius of 400 kpc; that is, 400 kpc marks the 
separation between the inner and outer power laws. 
Since all of our mass measurements are at radii smaller 
than this, our prof i le pro bes precisely the range where 
ISand. Treu fc EllisI ^200^ find a shallow inner profile. 

We believe this is due to the nat ure of their lens model. 
iBartelmann fc MeneghettH l)2004|) have shown that flat 
cores are indeed required for spherical and axisymmetric 
mass distributions, but if small deviations from axisym- 
metry are introduced (such as a mild ellipticity in the 
lensing potential, or a small external shear), these same 
systems are consistent with the cuspy profiles seen in 
simulations. 

3.3. Central density slopes and ACDM 

The confidence contours in Eigure El shown for the 
RCC subset, show a range of values for the central den- 
sity slope. This can be seen at the bottom of Eigure 0] 
we find that —1 < a < —2. This distribution ap- 
pears somewhat steeper than CDM predictions (the gray 
band). Using symmetrized errors, we find a weighted 
mean of (a)^ = -1.685 ± 0.077. The fact that this 
is marginally steeper than the CDM band is not sur- 
prising, since the contraction of cooling baryons is ex- 
pected to deepen the central cluster potential, and thus 
steepe n the central d ensity profile JHcnn awi & Ostrik eii 
1120021) : lGnedin et alJ^^2004^) : but see also lLoeb fc Peebles! 
l)20 03l)'). None of the clusters in the sample i s consistent 
with the flat or nearly-flat cores found by ISand etldl 
l|2004li) . We note, however, that only one cluster - 
MS2137 - appears in both samples. 

One can interpret these results in at least two ways. It 
is thought that baryonic processes could act to steepen 
the mass profile as energy is removed from the cluster 
core via relaxation and radiative cooling effects. The 
adiabatic contraction of the core baryons deepens the 
potential well, causing the d ark matter profile to steepen 
(jHcnnawi fc Ostrikeill200'^ . Another possibility is that 
merger tree hysteresis in the formation and evolution of 
the cluster, distinct from baryonic physics, dictates the 
current value of slope of the c ore dark matter profile 
()Ma fc de Bovlan-K olchin"200 ^ . The question of which 
influence - adiabatic contraction or accretion hysteresis 
- is dominant in setting the core profile, and the degree 
to which heating from a central AGN may play a role 
JOmma ct al. 2004; Omma fc Binney 2004), may remain 
unsettled until numerical experiments are able to self- 
consistcntly simulate the gas dynamics and cooling in 
conjunction with the gravitational evolution. 



4. CORE PROFILES AND SIDM 

We find that our mass profiles are adequately de- 
scribed by the theory of structure formation in an ACDM 
universe, so we can use them to bound the strength 
of the dark matter parti cle-particle self-interaction. 
iSpergel fc Steinhardl l)2000j) originally proposed SIDM 
as a solution to the core problem on galaxy scales. 
DSSW used numerical experiments to show that a self- 
interaction cross section ctdm in the range 0.5 to 5 cm^ 
g^^ (~ 1 to 10 b(GeV/c^)~^) provided the requisite flat- 
tening for galaxy-sized dark matter haloes. As they have 
noted, a cross section of this size is probably inconsis- 
tent with the distribution of halo shapes in cosmological 
simulations (see also fYos hida et al, (2000)). They also 
point out the difficulty inherent in core profile compar- 
isons between simulations and observations, noting the 
i ndeterminate effect of baryon p hysic s . 

lArabadiis. Bautz fc Garmirel (j2002) showed that the 
mass profile for the relaxed cluster MS1358 was inconsis- 
tent with ctdm > 0.1 cm^ g~^. This w as based on a com- 
parison with the lYoshida et all 1)20001) simulations which 
provided an informal relationship bet ween cluster core 
size and self-interaction cross section. 'Marke vitch et all 
(■2003.) noted the need for a sample of clusters to place 
a robust limit on the cross section, and derived a less 
restrictive limit of 1 cm^ g~^ based on a small offest be- 
tween the weak lensing and X-ray peaks of a supersonic 
subcluster in IE 0657-56. 

In order to gauge the dark matter self-interaction 
strength, we determine the maximum core size of each 
cluster in the sample. Although we see no evidence for 
a constant density core in any cluster, we calculate the 
size of such an unresolved core that would be consistent 
with the data. We m odel each core as a softened isother - 
mal profile sphere t Arabadiis. Bautz fc Garmirel l200'^ . 
fitting for its size and density normalization: 



p{r) 



2-KG{r^+rl) 



(4) 



where a is the (constant) ID velocity dispersion and 
is the core radius. Since our profiles are determined in 
cumulative mass, rather than density, we fit for Tc and 
Mg using Mr = Mg {x — tan^-'^ x), where x — r/r^ and 
the mass scale Mg = ^a^rdG. 

To limit the variance of the fitted parameters which 
are due to a non-hydrostatic core, we do the analysis 
only for the RCC subset (see Table Confidence con- 
tours in core radius and mass scale are shown in Fig- 
ure 13 A similar relationship obtains between core ra- 
dius and two-phase plasma to that found for core slope: 
two-phase emission models result in smaller core sizes or 
upper limits. As in the case of core slope, the effect is 
most pronounced in A2029. Thus the limits derived from 
Uniphase core models can be considered to be conserva- 
tive estimates. 

The core radii of the RCC subset are shown at the 
top in Eigure 0] (The corresponding core slopes for the 
same clusters are shown below.) The dashed lines show 
core sizes expected for three values of the dark m atter 
self-interaction cross section ifYoshida et al.ll2000l) . The 
DSSW cross section required to solve the core prob- 
lem - 0.5 to 5 cm^ g"-'^ - is clearly ruled out. Us- 
ing symmetrized errors, the weighted mean core size is 



(re) = 35.9 ± 15.4 kpc, corresponding to self-interaction 
cross section ctdm ^ 0.1 cm^ g~^. We place a 99% upper 
limit of 65 kpc on the cluster core size, corresponding to 
just under 0.1 cm^ g~^. 

A few caveats must be mentioned here. The first is 
that this upper limit is derived from a calibration with 
simulations which did not include the effects of baryon 
physics in the core evolution. Presumably the adiabatic 
contraction of the core from radiative losses could coun- 
teract some effects of dark matter self-interaction, and 
could be consistent with the observations, provided that 
the contraction timescale is comparable to the mean 
time between dark matter particle scatterings (to pre- 
vent cooling- induced core collapse), and that the self- 
interaction is not so vigorous that gravitational core col- 
lapse proceeds in a Hubble time. However it is not clear 
that the baryons only act to hasten core contraction. As- 
trophysical processes associated with AGN or star forma- 
tion, such as jetc/ICM interactions, supernovae and stel- 
lar winds, may be required to match observablcs such as 
the "excess entropy" in groups and clusters Frenk ( 200^, 
although the issue of feedback is beyond the scope of this 
paper. 

A second point is that the sel f-interaction cross se ctions 
in the calibration simulations ijYoshida et al.ll2'00GD were 
assumed to be velocity independent. DSSW pointed out 
that a cross section which scaled as might ameliorate 
the discrepancies between SIDM cluster simulations and 
observations. A functional dependence such as aoM — 
aQ{v / va)~^ might reduce the effects of self- interaction on 
cluster scales, although early results are not promising 
(N. Yoshida, private communication). 

Finally, if the properties of t he core are in fact governed 
by clu ster accretion hysteresis (jMa fc de Bovlan-Kolchiiil 
120041) , measurements of cluster core properties will prob- 
ably tell us little until we have a sufficiently large obser- 
vational sample that we can explore the relevant physics 
through statistical analyses. 

5. SUMMARY 

We have presented an analysis of 10 galaxy clusters 
observed with the Chandra X-ray Observatory. We have 
carefully modelled their cores, utilizing a second plasma 
component where required. We have derived mass pro- 
files for each of them, and for the RCC subset we find 
that ACDM provides an adequate description of their 
cores. We have placed upper limits on the size of any 
constant density core in each cluster, and derive a 99% 
upper limit of just under 0.1 cm^ g^-'^ — 0.2 b(GeV/c^)~^ 
for the dark matter self-interaction cross section. 
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Fig. 1. — Mass profiles (top) and core power law slopes (bottom) of each cluster in the sample. We show profiles for Uniphase (left) and 
two-phase (right) models of the core plasma (except for Hydra A; see text). The 'preferred' label appears on the model with the higher 
MCMC significance. 
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Fig. 2. — Confidence contours of power law fits to the cores of the RCC subset. The left/right plot in each pair is for uniphase/two-phase 
models of the core plasma. The model with the higher MCMC F-test significance is labeled 'preferred'. 
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Fig. 3. — Confidence contours of softened isothermal sphere fits to the cores of the RCC subset. The left/right plot in each pair is for 
uniphase/two-phase models of the core plasma. The model with the higher MCMC F-test significance is labeled 'preferred'. 
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Fig. 4. — Summary of the softened isothermal sphere core radii (top) and power law slopes (bottom) for fits to the RC C subset. The 
relatio n between the self-interaction cross section of dark matter particles aoM and the cluster core radius Tc, taken from lYoshida et alj 
120001) . is d enoted by d ashed lines. The inner slope predictions of numerical CDM exp eriments is shown by a gray band bracketed by the 
values from INavarTolT'r cnk & White 1 1996, 1997) at a = -1 and lMoore et alj 11999bl) at a = -1.5. 
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TABLE 1 
Galaxy Cluster Sample 



cluster 
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core plasma 
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RCC? 
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0.0522 


Uniphase 




no 


Abell 1795 
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Uniphase 
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yes 
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Uniphase 


0.483 


yes 


ZwCl 3146 


0.2906 


two-phase 


0.989 


yes 
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